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Abstract. The developments in this paper are concerned with nonholonomic 
field theories in the presence of symmetries. Having previously treated the case 
of vertical symmetries, we now deal with the case where the symmetry action can 
also have a horizontal component. As a first step in this direction, we derive a 
new and convenient form of the field equations of a nonholonomic field theory. 
Nonholonomic symmetries are then introduced as symmetry generators whose 
virtual work is zero along the constraint submanifold, and we show that for 
every such symmetry, there exists a so-called momentum equation, describing the 
evolution of the associated component of the momentum map. Keeping up with 
the underlying geometric philosophy, a small modification of the derivation of the 
momentum lemma allows us to treat also generalized nonholonomic symmetries, 
which are vector fields along a projection. Such symmetries arise for example in 
practical examples of nonholonomic field theories such as the Cosserat rod, for 
which we recover both energy conservation (a previously known result), as well 
as a modified conservation law associated with spatial translations. 
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1. Introduction 

Knowledge of symmetry is fundamental to the understanding of mechanical systems 
and field theories, whether classical or not. Results such as Noether's theorem, the 
Marsden-Weinstein reduction procedure, etc. all attest to this. In the presence of 
nonholonomic constraints, however, the concept of symmetry becomes somewhat more 
subtle: Noether's theorem (see [10] and the references therein) for example no longer 
holds automatically. 

Surprisingly, there is a more sophisticated analogue of Noether's theorem, going by the 
name of the nonholonomic momentum equation [5,7], which plays a fundamental role 
in the reduction procedure for general nonholonomic systems. This equation describes 
the evolution of the components of the momentum map under the nonholonomic flow 
and its derivation relies heavily on the use of differential geometry for the formulation 
of the nonholonomic problem. 

In the case of classical field theories, a similar result was derived in [22] for a restricted 
class of symmetries, namely those symmetries whose infinitesimal vector field is 
vertical, or alternatively, symmetries which act trivially on the space of independent 
variables of the theory. Nevertheless, many physically interesting symmetries are 
not vertical, and an extension of this result to a more general class of symmetries is 
therefore needed: this is the subject of the present paper. Examples of non-vertical 
symmetries include translation in time, yielding conservation of energy, without doubt 
the prime example of a conserved quantity. 

Outline of this paper 

In this paper, we derive a generalized form of the nonholonomic momentum equation 
for two classes of symmetries. The first consists of symmetries which have a nontrivial 
horizontal component, while for the second class we restrict our attention to vertical 
symmetries whose infinitesimal generator is a generalized vector field, i.e. a vector 
field along the projection tt^o- More information on the application of generalized 
symmetries to differential equations can be found in [19]. 

In section |3] we derive a new form of the field equations for a classical field theory with 
nonholonomic constraints, which is particularly well-suited to the derivation of the 
momentum equation in section!?] These equations are equivalent to the nonholonomic 
field equations derived, for example, in [3,15,21] but they do not involve the Lagrange 
multipliers, and arc therefore closer in spirit to Holder's equations (sec [18]). We 
derive the nonholonomic momentum equation for two distinct classes of symmetries: 
firstly, for the case of a symmetry group whose infinitesimal generators arc protectable 
vector fields on the total space, and secondly, for a symmetry group acting vertically 
on the total space, but whose associated nonholonomic symmetries are generalized 
vector fields. 

This theory is applied in section to a number of examples. In section 15.11 we 
obtain a new class of nonholonomic symmetries for mechanical systems with nonlinear 
constraints; Benenti's nonholonomic system (see [2]) is chosen as a straightforward but 
relevant example and a number of new conservation laws are derived. In section [5T2"l we 
derive nonholonomic conservation laws for the nonholonomic Cosscrat rod, an example 
of a nonholonomic field theory (see [23]). These conservation laws are associated to 
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translations in time and spatial translations, respectively. The former is a horizontal 
symmetry, while the latter is described by a (vertical) generalized vector field. 

Finally, in the appendix, we elaborate somewhat further on the choice of reaction 
forces for a nonholonomic field theory. As the material in this section is not crucial to 
the developments in the remainder of the paper, we have delegated it to the appendix. 

Notation 

Throughout this paper, all geometric objects are assumed to be smooth. We will 
denote the contraction of a vector field V and a form a both as iyct or as VJ a, and 
we will use these two notations interchangeably. The Lie derivative is denoted by „£? . 
We will frequently use the notation /\ k (M), where M is a manifold, to denote the 
fc-fold exterior product of TM with itself. Thus, /\ k (M) is a bundle over M whose 
sections are fc-forms on M. 

2. Classical field theories 

As is customary in the geometric treatment of classical field theories (see [8, 10, 13,20] 
and the references therein) , fields arc modelled as sections of a fibre bundle it : Y — > X 
of rank m, where the base space X is an (n + l)-dimcnsional oriented manifold (with 
volume form rj), and the total space Y has dimension n + m + 1. Local coordinates 
on X are denoted by (x^), \i = 0, . . . , n and are supposed to be such that the volume 
form rj can be locally written as 

n = d n+1 x := dx Q A • • • A dx n . 

In addition, we assume a local system of bundle coordinates (x^ 7 y a ), fi = 0, ...,n, 
a = 1, . . . , to, on Y to be given. 

Over Y, there exists a tower of jet bundles 
■ ■ • — > J A n > JV — ► K 

The elements of J k ir are equivalence classes of local sections of it, where two sections 
are said to be equivalent at a point x of X if their fcth-order Taylor expansions agree 
at that point. We denote the equivalence class of a section as j£<fi. More information 
on jet bundles can be found in [20]. We will mostly only need the first-order jet bundle 
J 1 n; the second-order jet bundle will make a brief appearance in section [5T2l (see for 
instance [12] for a review of mechanics on higher order jet bundles). 

We will need as basic geometric tools the concepts of contact forms and jet 
prolongation of vector fields. A contact m-form on J k ir will be any TO-form 9 satisfying 
(j <f>)*0 for every local section <j> of it. The set of contact 1-forms on J lr K defines 
a distribution ©cartam called the Cartan distribution. Using this distribution, it 
is easy to define the 1-jet prolongation of a vector field £y on Y as the unique 
vector field j £y projectable on £y and preserving the Cartan distribution, i.e. 
^fj^ Y ©Cartan Q ©Cartan- In coordinates, if (, Y = C M (a;,2/)gfjr + £, a ( x , y) gSr> then 
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If we adopt the coordinate systems on X and Y as above, then J 1 it is equipped 
with an induced coordinate system y a ;y^), /U = 0, . . . ,n, a — 1, . . . , m. With 
respect to this coordinate system, the projection tti : J 1 ^ — > Y is given by 
7r 1)0 (^,y«;y°) = (a^.y ). 

For later use, we mention the existence of a distinguished vector- valued (n + l)-form 
S v on J 1 ^, which is called the vertical endomorphism. In coordinates, S v reads 

S„ = (dy a -y»d^)Ad n z M ®^, (1) 



where 



For the purpose of this paper, a Lagrangian will be a function L on J 1 it. Given a 
Lagrangian L, we define the Poincare-Cartan (n + l)-form as 



dL 



6 L := S;(dL) + L V = — — (dy a - y'dc") A d"^ + Ld 



n+l 



,r. 



We also put fix := — dOj,, to which we refer as the Poincare-Cartan (n + 2)-form. 
For more information on the geometry of jet bundles, see [20]. 



The Euler- Lagrange equations 



Let there be given a Lagrangian L. The dynamics of the field theory is described 
by the Euler-Lagrangc equations associated to L; they express that the field is an 
extremum of the following action functional: 

5(0)=^L^,0 a (x),^d" +1 x. (2) 

By varying the action S with respect to a variation V, we obtain after integrating by 
parts 

"-sW'^L-X^-s^)™* 1 * (3) 

where $ e is a finite variation associated to V. Here, a variation of a field <f> over an 
open set U C X is a 7r-vertical vector field V defined on an open neighborhood of 
4>(U). The associated finite variation <!> e is nothing but the flow of V. Note that the 
composition $ e o <p is again a local section of tt. 

In coordinates, a section of 7r is an extremum of @ if and only if it satisfies the 
familiar Euler-Lagrangc equations, given by 

There exist various intrinsic formulations of these equations (see [11] for an overview), 
of which wc mention just one. It can be shown by a straightforward coordinate 
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calculation that the Eulcr-Lagrangc equations are equivalent to the following set of 
intrinsic equations: 

C/V)*(% n i)=0 for all W £ X(J 1 tt). 

In what follows, we will be mostly interested in similar equations for field theories 
which are subject to nonholonomic constraints. This is the subject of the next section. 

3. The nonholonomic field equations 

A nonholonomic field theory is given by the specification of three objects (see also [21]): 

(i) a Lagrangian L : J lr K — > R; 

(ii) a constraint submanifold C e — > J 1 ^, such that the restriction of the projection 
( 7r i,o)|c defines a sub-bundle of nifi : J it — > Y; 

(iii) a bundle of reaction forces F, where the elements 4> of F are (n + l)-forms defined 
along C, i.e. maps from C C J 1 n to (J 1 ^). The elements of F have to satisfy 
the following requirements: 

(a) $ is n- horizontal, i.e. $ vanishes when contracted with any two 7Ti-vertical 
vector fields; 

(b) $ is Tcontact, i.e. (j 1 </>)*<I> = for any section <fi of n. 

It can be shown that any element <& of F is of the following form: 

$ = ^(dj/ -2/«d^)Ad"^, (4) 
where A 1 ^ are functions on J lr K. 

For the sake of simplicity, we will assume C to be defined by the vanishing of k 
functionally independent functions ip a on J lr K. Furthermore, we will assume that F 
is globally generated by I generators $ K of the following form: 

^ = A^(dy a - ytdx-) A d"^ (« = 1, . . . , I). 

In practice, the dimension / of F will be equal to the codimension k of C. There seems 
to be no a priori reason for supposing that k = I. In most cases, however, F will be 
determined by C through application of the Chctaev principle, described in remark 
below. 

In the nonholonomic treatment of constraints, a special role is played by infinitesimal 
variations of the fields which are compatible with the constraint, as in the following 
definition. 

Definition 3.1. A variation V of a field if) (taking values inC, i.e. such that j 1 ^ ClC) 
defined over an open subset U with compact closure is admissible if 

(jWO' 1 ^ *) = for all® g F. (5) 

By varying the action S with respect to an admissible variation V, we obtain after 
integrating by parts again If the variations V were arbitrary, then this would 
immediately yield the Euler-Lagrangc equations. However, this is not the case as 
the variations have to be admissible. Additional reaction forces will therefore appear 
in the Euler-Lagrangc equations, whose role it is to constrain the solution <\> to the 
constraint submanifold. 
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Definition 3.2. A local section <f> of -k , defined on an open subset U C X with compact 
closure, is a solution of the nonholonomic problem determined by L, C, and F if 
j <p(U) C C and (0) holds for all admissible variations V of <f>. 

It follows from ([3]) that a local section </> is a solution of the nonholonomic problem if 
it satisfies the nonholonomic Euler- Lagrange equations: 



dL d dL 



dy a dx^ dy® 



O'V) = Aa«^r O'V) and ^(jV) = 0. (6) 



Here, A QK are unknown Lagrange multipliers, to be determined from the constraints. 
An intrinsic form of these equations is derived below in theorem 13.51 but first we need 
the following technical results. 

Lemma 3.3 (lemma 3.2 in [13]). Let W be a vector field on J l ix. If (p is a section 
of 7r and if either W is tangent to the image of j 1 ^ or if W is %ifl- vertical, then 
(j 1 4>)*(iwn L )=0- 

Now, let be a section such that the image of j 1 ^ is a subset of C and consider a 
vector field W which is tangent to the image of j 1 ^, i.e. there exists a vector field 
w on X such that T x j 1 4>(w(x)) = W(j^.(f)) for all x G X. One can follow a similar 
reasoning as in the proof of lemma 3.2 in [13] to show that 

(jV)*(wj<t>) = ™j((jV)*$) 

for any $eF. Since <!> is 1-contact, the right-hand side of this expression vanishes. On 
the other hand, if W is tt\ ^-vertical, it follows automatically that (j 1 </>)*(M / J $) = 0. 
We have therefore proved the following lemma: 

Lemma 3.4. Let <p be a section of ir such that j^<j> £ C for all x £ U C X. If either 
W is tangent to the image of j l 4> or W is iri^- vertical, then (j 1 0)*(H / J <!>) = for all 

$ e F. 

Henceforth, we shall call any vector field W on J 1i k admissible with respect to a section 
cj) of 7T if (>'V)*(WJ $) = for all $eF. 

Theorem 3.5. Let <f> be a section of it. Iflmj 1 ^ C C, then the following assertions 
are equivalent: 

(a) (p is a stationary point of the action (0) under admissible variations; 

(b) <j) satisfies the Euler-Lagrange equations §Bj); 

(c) for all vector fields W on J 1 it such that (j 1 (j))*(W} <&) = for all $ E F, 

(j 1 4>y(WASl L ) = 0. (7) 

Proof: Let us first prove the equivalence of (jlj) and (jg). For arbitrary, not necessarily 
admissible variations, the following result holds (this is equation 3C.5 in [13]): 

^-s(^) =- f ( J 1 cj>nj 1 v_\n L ). 



de 



e=0 



U 



For admissible variations, from hypothesis (a), we have 



V 
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Now, we may multiply V by an arbitrary function on X and this result will still hold 
true. The fundamental lemma of the calculus of variations therefore shows that 

V 1 <t>)*V 1 vm L ) = o, (8) 

for all admissible variations V defined over U. By using a partition of unity as 
in [13], it can then be shown that ((SJ) holds for all n- vertical vector fields V such 
that (j 1 V)A $ = for all $eF. This expression is equivalent to 0: to see this, take 
an arbitrary vector field W on J r ir such that (j 1 (j))*(WA $) = for all $ e F. The 
vector field W can be decomposed as the following sum (to be considered along the 
image of j V) : 

W = w\\ +j 1 V + v nit0 , 

where w\\ is tangent to the image of j (f>, j V is the prolongation of a 7r-vertical vector 
field V, and v ni is a Tr^o-vertical vector field. Using lemma [3"^fl we have that 

ti 1 <i>y(j 1 vi$) = (j i <i>nwi$) = o, 

and from lcmmaGU we get (jV)*(^J Hi) = {j 1 <t>)* tfV 1 Q.l)- The right-hand side 
of this equation vanishes since j x V is admissible, and therefore we conclude that W 
satisfies (J7J. 

The equivalence of jbj and Jcj is just a matter of writing out the definitions. In 
coordinates, the left-hand side of Q reads (for a prolongation of a vertical vector field 
V) 



and this holds for all admissible variations V. Therefore, if satisfies ([7]), then there 
exist functions X aK such that 



d dL 



9y a dx^ 1 

The converse is similar. □ 

We see from the proof of this theorem that only vertical vector fields yield nontrivial 
results for Q. 

Remark 3.6. The bundle of reaction forces F is commonly derived from the constraint 
submanifold C through application of the Chetaev principle (see [18], as well as [21] for 
an extension to the case of field theories). If the constraint submanifold is given as the 
zero level set of functions cp a , then according to this principle, F is locally generated 
by the following forms: 

$ tt := 5*(d^) = §^(dy a - ytdxn A d»V 

In the past, there has appeared some criticism over the use of the Chetaev principle 
(see [16]), and as we shall see in the appendix, for classical field theories the Chetaev 
principle sometimes has to be modified. 
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4. The nonholonomic momentum equation 

In this section, we derive the nonholonomic momentum equation, the nonholonomic 
counterpart to the well-known theorem of Noether. More in detail, we prove that 
for every nonholonomic symmetry there exists a certain partial differential equation, 
which reduces to a conservation law when the constraints are absent. In proving the 
momentum lemma, many different starting assumptions can be made, and we study 
two different setups: 

(i) in section 14.21 we assume that the symmetry group G acts by bundle 
automorphisms on Y . In particular, we allow for the fact that G acts nontrivially 
on X as well. In this case, nonholonomic symmetries are projectable vector fields 
on Y. 

(ii) in section 14. 3[ we model nonholonomic symmetries as vector fields along the 
projection ni t o- For the sake of simplicity, we assume in this case that the 
action of G on Y is vertical, i.e. G acts trivially on X. This setup is therefore 
complementary to the one described before, but this case can probably be 
extended even further. 

4-.1. The nonholonomic momentum map 

Let G be a Lie group acting on Y by bundle automorphisms; i.e. there exist smooth 
actions fiGxY^F and $:GxI^I such that 7r(I>(#, y)) = J>(g, n(y)) for all 
g G G and y £ Y. We use the following abbreviations: Q g := <&(<?, •) and $ g := ■)■ 
The Lie group G acts on J 1i k by prolonged bundle automorphisms, i.e. 

.,' <!>,(./;,,; --l!V : - 

Let us now assume that G leaves invariant L, C, and F: 

3 1 & g (Lq) = bn, j^CjcC and (j l <& g )* F C F 

for all g £ G. 

We consider first the vector bundle q f over Y, defined as follows. Denote by g F (y) 
the linear subspace of g consisting of all £ £ g such that 

i 1 er(7)Ji ? = for all 7 eCn7rrj(y), (9) 

where £y is the infinitesimal generator of the action corresponding to £, that is, 

&(y) = T7^exp(tf)(2/) 

at t=o 

We assume that the disjoint union of all Q F (y), for all y £ Y can be given the structure 
of a vector bundle q f over Y. 

To any section £ of g F , one can associate a vector field £y on Y according to the 
following prescription: 

fr(») = [e(v)] y (»). (io) 
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Definition 4.1. A nonholonomic symmetry is a section £ of g F such that the 
associated vector field £y is ir-projectable; i.e. there exists a vector field £x on X 
such that Tir o £y = £,x ° 7T - 

Definition 4.2. A horizontal nonholonomic symmetry is a section £ of g F which is 
constant, i.e. = t;(y') for all y, y' G Y. 

We may identify horizontal nonholonomic symmetries with elements £ G q such that 
£ = Q F (y) for all y G Y, that is, 

j^y (7) J F = for all 7 G C n tt^ (y) and for all y G Y, 

We now define the nonholonomic momentum map as the map J ' : C — > /\ n ( J 1 7r) <E) 
g F , constructed as follows. Let £ be any section of g F (for the construction of the 
momentum map it does not matter whether the associated vector field is projectable 
or not) and put 

where £ is the vector field associated to £ according to (|T0|) . Given this definition of 
, we define J n h - by the following rule: 

(J n -\e) - Jf- h \ 

If £ is a nonholonomic symmetry, then the prolongation j £y of the associated vector 
field is admissible. This is proved below in corollary 14.41 

Lemma 4.3. Let £ be a section of g F . For a fixed y £Y, put £ := and consider 
any 7 G 7rf q(^) ^ ^. Then there exists a ni^-vertical vector t> 7 G T 7 J 1 ^ swc/i that 

/M7) = j 1 ?F(7) + w 7 - 

Proof: Recall that j^y is the prolongation of the fundamental vector field £y 
associated to the (fixed) Lie algebra element £. The lemma follows from the fact 
that T 7 7r l! o(j 1 |y( 7 )) = £ Y (y). On the other hand, r 7 7ri, (j 1 fr(7)) = My): but b Y 
definition, |y(y) = £y (y). □ 

Corollary 4.4. Let ^ be a nonholonomic symmetry. Then the prolonged vector field 
j £y is admissible with respect to any section <f> of it. 

Proof: Take 7 G C, put y = 71-1,0(7), an d let £ = £(?/)■ From definition 14. 11 we gather 
that for all $ 6 f, j^y (7)J $(7) = 0. But j 1 ^ (7) is equal to j^y (7) up to a vertical 
vector, while <f> G F is semi-basic. Therefore, 

J^y( 7 )J $(7) = 0. 

As this holds for all 7 G C, we conclude that j t~y is admissible with respect to any 
section. □ 

Remark 4.5. In the case of mechanical systems with linear constraints, the bundle 
g F , defined as above, coincides with the bundle S used in [5]. o 
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4-2. The nonholonomic momentum equation 

Theorem 4.6. If <f> is a solution of the nonholonomic field equations Q), then for 
any nonholonomic symmetry £ the associated component of the momentum map J^ h ' 
satisfies the following nonholonomic momentum equation: 

(jV)*(djf- h -) = yvr (^ iY (L V )). 

Proof: Let f be a nonholonomic symmetry. Recall that the prolongation j 1 ^ of the 
associated vector field is admissible with respect to any section of it. Therefore, since 

pulling back along a solution of the nonholonomic Euler-Lagrange equations gives us 
(jV)*(dJf h ') = (iV)* (^V©l). (12) 

For the sake of notational convenience, let us say that two forms a and /3 on J 1 ^ 
are equivalent (denoted by a ~ /3) if they agree up to a contact form, i.e. a ~ (3 iff 
a = (3 + 9, where 8 is contact. This is equivalent to saying that (j 1 <f>)*a = (j 1 ^ 1 )*^ for 
all sections <fi of ir. We then have that 

^£ y ®l = Se^y {S;(dL)) + X. Hy (Lr,) . 

In coordinates, the first term on the right-hand side becomes 

s> jlir (s;(dL)) = se Hy {^Mt - tM) a d"* M ) 

BT 

^^ V(V " J,X)A<1% 
BT Bf a 

= ^^r(dy b - y h Ax v ) A d"^ ~ 
By* <V y y " ' ^ 

and hence vanishes when pulled back along a prolongation of a section. Here, we've 
used the fact that the prolongation j 1 ^ can locally be written as follows: 

and that the Lie derivative of a contact form with respect to such a vector field is 
again a contact form. Note especially that the coefficient £^(x) does not depend on 
y: this is a consequence of the projectability condition in definition 14.11 □ 

Corollary 4.7. If <f> is a solution of the nonholonomic field equations then for 
any horizontal nonholonomic symmetry £ the momentum map J^ l h ' is a conservation 
law for the nonholonomic dynamics, that is, 

(jV)*(dJ?- h -) = 0. 
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4-3. Vertical vector fields along the projection 

Not all nonholonomic symmetries can be modelled as sections of g F ; as we shall see, an 
important class consists of sections of the pullback bundle irl g F , i.e. maps £ : C — > g F 
such that £(7) £ 9 F (tti.o(7)) for all 7 € C Such a section can be represented as 

£=t A (x>*,y a ,i£)e A , (13) 

where £ A (x li ,y a ,y?i) are locally defined functions on C and {e^} is a basis of g. 

For this treatment, we will make one important simplification: we assume that the 
symmetry group G acts vertically. It should noted, though, that this is probably not 
a fundamental restriction: it is likely that this approach can be generalized further. 

Our treatment is therefore somewhat restricted compared to the previous paragraph. 
On the other hand, by allowing as in (| 13[) sections whose coefficients are functions on 
C rather than just on Y, we are able to describe symmetries which do not fit into the 
framework of the previous paragraph. Examples will be given in section [5] 

Definition 4.8. A generalized nonholonomic symmetry is a section £ of the pull-back 
bundle tt\ g F . 

A generalized nonholonomic symmetry £ induces a vertical vector field £ along 7Ti o s 
defined as follows: for all 7 £ C. 

1(7) = [£(7)] Y (V) € T v Y > whcre V = ^oW- 

There exists a generalized notion of prolongation for vector fields along 71-^0 (see [20]). 
In coordinates, if £ = £ a (x^, y b , y^)-^, then 



5 dy a dx» dy^ W 01 C dxn dx» dy bV » dy b v V ^' 



Note that j -1 £y is a vector field along 7r 2 ,i, as its coefficients depend on the second-order 
derivatives y^. 

Having defined the prolongation j'^y as above, we now define the associated 
component of the momentum map as 

This definition is formally identical to (|11|) . but strictly speaking, the momentum map 
is now an n-form along the projection 1*2,1- By pull-back, such a form induces an 
n-form on J 2 ir. Similarly, the concept of Lie derivation can be extended in a natural 
fashion to the case of vector fields along the projection (see [20]). 

Theorem 4.9. If <f> is a solution of the nonholonomic field equations ([?]), then for any 
generalized nonholonomic symmetry £ the associated component of the momentum map 
Jf satisfies the following nonholonomic momentum equation: 

O' 2 0)*(dJf- h -) = (j^y (^H Y (LV)) • (14) 
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Proof: The proof is similar in spirit to that of theorem 14.61 but there are some 
additional technical difficulties that need to be taken into account. For any solution 
<fi of the nonholonomic field equations, we have 

(iV)*(dJ|' h -) = (j 2 d>n^® L ) + (fcf>y(i jliY n L ). (15) 

For the first term on the right-hand side, one can check that 

^v( d y a - y a ^) - ^( d y b - y>") + - y>^ 

a contact form on J 2 tt. From this, it follows that the Lie derivative of Ql is given by 

The second term on the right-hand side of (|15p is zero. This is essentially a consequence 
of the fact that <f> is a solution of the field equations and can be proved using the 
following observation. Take x S X and consider a vector field W on J 1 ^ such that 
W(jl(j>) = fiyijl^)- Then W(jl<j>)}$ = for all $ e F, and without loss of 
generality, we may choose W to be such that this equality holds in a neighbourhood 
of jx</>. We now have that 

C?V)*(Vfrfi£)(aO = (j X <t>y{iw&L){x). 

As W is admissible with respect to <fi, we conclude that the right-hand side is zero, 
and this in turn implies the vanishing of the second term in (|15j) . □ 

5. Examples 

5.1. Mechanical systems with nonlinear constraints 

As an illustrative example, we take a variation of Benenti's system (see [9, 14] and the 
references therein), which describes two point masses moving on a horizontal plane 
whose velocities arc constrained to be parallel. The configuration space of this system 
is Q := R 2 x K 2 , and so we put X = E, while Y = K x Q. The projection n is the 
projection onto the first factor. The Lagrangian for this system is 

L = j{±l + yi) + j{x 2 2 + y 2 2 ). 

Strictly speaking, L is a function on TQ, but since J 1 it is isomorphic tolx TQ, we 
view L as a (time- independent) Lagrangian on J 1 ^. The constraint is 

ip = XxVl ~ X2Vl = 0. 

This is a nonlinear constraint, and determines a submanifold C of J 1 ^. If we assume 
Chetaev's principle to hold, the bundle of reaction forces F along C is generated by 
the following one-form: 

$ = 2/20xi + iiO y2 - yiO X2 - x 2 Vl , 
where the contact one-forms 9 q a are defined as 9 q a = dq a — q a dt. 
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(a, (3, 7, 5) 




(1, 1, 0, 0) 


±ixi + yiy\ 


= 


(0, 0, 1, 1) 


X2X2 + wh 


= 


(1,0, 0,-1) 


x\ii - y 2 y2 


= 



Table 1. Conservation laws for the Benenti system. 



Consider now the obvious action of the Abelian group R 4 on Y by translations: 
(01,61,02,62) • (t\x\, y\,X2,V2) = (t;xi + a\,yi + b l7 x 2 + a 2 ,y 2 + h)- This action 
fulfills all the necessary conditions needed for the momentum lemma: it is vertical, 
and leaves invariant the Lagrangian L. the constraint submanifold C, and the bundle 
of reaction forces F. 

It is easy to check that the contraction of a vector field £ of the following form 

d d d d 

£ = axi- \-/3yi- r-7±2~ h Sy 2 Tr (a, (3, 7, S G R) 

0x1 oyi 0x2 ay 2 

with $ will vanish along C if a+5 = (3+j. In this case, £ is a generalized nonholonomic 
symmetry of the kind defined in definition 14.81 (Strictly speaking, the term 
"generalized nonholonomic symmetry" refers to the section £ := (axi, (3yi, 7x2, Si/2) 
of^ of /.) 

For the sake of convenience, let us take a = j3 = 1 and 7 = 5 = 0. The other cases 
are similar and the results are summarized in table [T] In this case, the prolongation 
j £ of £ is the vector field along 7^1 given in coordinates by 

1 r . d . d d .. d 

ox\ oyi 0x1 oyi 

The component of the nonholonomic momentum map then becomes 

Jf h -=j 1 OQL = m{x\ + yl), 

while the right-hand side of the momentum equation (| 14[) is Jzf^^Ldt) = m(iixi + 
yij)i)dt. The momentum equation hence reduces to iiXi + yijji = 0, which can also 
be verified using the equations of motion. 



5.2. The nonholonomic Cosserat rod 

The nonholonomic Cosserat rod is an example of a nonholonomic field theory studied 
in [23]. It describes the motion of a rod which is constrained to roll without sliding 
on a horizontal surface. This theory can be studied using the bundle tt : Y — > X, 
where X=[0,/]xl (space and time) and V^IxtfxS 1 . with bundle coordinates 
(s,t;x,y,8). 

Its Lagrangian is given by 
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with k 2 = (x") 2 + {y") 2 , while the constraints are given by 

x + My' = and y - R9x' = 0. (16) 

Here, p, a, j3, K, and R are real parameters. The field equations associated to this 
Lagrangian are given by 

px + Kx"" = A 

py + Ky"" = fx (17) 
ad - (39" = R(\y' - fix'), 

where A and p are Lagrange multipliers associated with the nonholonomic constraints. 
These equations are to be supplemented by the constraint equations (fPS)) . 

For future reference, we note that the bundle of reaction forces is generated in this 
case by the following forms: 



= (dx-xdt) A ds + Ry' (d9 - 9dt) Ads 
$ 2 = (dy — ydt) A ds — Rx'(d9 — 6dt) A ds. 



In the absence of nonholonomic constraints, this model is subject to the usual 
symmetry actions such as translations in time, global translations, and global rotations 
(see [23]). As we shall now show, some of these persist in the nonholonomic case. 

Observe that the Lagrangian is of second order. At the cost of sacrificing physical 
relevance, one may also put K = to obtain a purely first-order theory. We will 
not do so here, as the derivation of a nonholonomic momentum lemma for a second- 
order field theory proceeds exactly as above, up to a few minor modifications. The 
nonholonomic momentum map J n h ' is now defined on J 3 7r, and the nonholonomic 
momentum equation hence becomes 

CfV)*(dJf' h -) = (j 2 <t>TJ? jH (Lr,), (19) 

where is a solution of the nonholonomic Euler-Lagrange equations. This equation 
holds both for "genuine" symmetries as well as for symmetries along the projection. 



5.2.1. Translations in time Consider the action of R on X by translations in time 
defined by the map $ : K x X — > X, with <f>(a, (s, t)) = (s,t + a). As the bundle n 
is trivial, this action naturally induces an action on Y , and, by prolongation, also on 
J lr K. Clearly, this action is not vertical; the fundamental vector field associated to a 
Lie algebra clement £ G R is given by 

i^=4 (20) 

This vector field is a nonholonomic symmetry: it is 7Ti-related to a vector field on X 
and its contraction with the elements of F vanishes along C: for the n-form $ x defined 
in (fTS)) . we have 

-^J^ 1 = -(x + R6y')ds, 
which vanishes on C, and a similar result holds for $ 2 . 
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Note that the nonholonomic symmetry (|20p is a horizontal symmetry (definition [ 
in general, this will not be the case. In this special case, however, we have that 
Jtfji£ Y (Lr]) = 0, expressing the infinitesimal invariance of the Lagrangian. The 
momentum map now becomes 

(f(f>)*J^ h - = \-Kx"'x + -Ky"'y + 09 + K(x"x' + y"y')] dt + Eds, 
where we have introduced the energy density 



£ = f (i 2 + v 2 ) + f^ 2 + 1((^") 2 + {y") 2 ) + §(*T- 

The nonholonomic momentum lemma then states that d(j 3 </>)* Jf' 11 ' = 0. This 
equation expresses local conservation of energy; by integrating over a hypersurface 
of constant s, we may then obtain a law expressing global conservation of energy. 

Remark 5.1. This case was also treated in [23] by use of a different method, to which 
we refer for further details. Note however, that the method used in that paper is less 
general as it is only valid for horizontal nonholonomic symmetries. o 



5.2.2. Spatial translations Consider the action of R 2 x S 1 on Y by translations; i.e. 
for each (a, b, ip) we consider the map $( a ,6, v ) : (s, t, x, y, 9) i— > (s, t, x + a, y + 6, 9 + tp). 
Let £ = (vi,V2,vg) be an element of the Lie algebra of I 2 x S 1 . The corresponding 
fundamental vector field is given by 

d d d 

CY =Vl-£- + Vi— + V0-. 

ox ay ov 

When no constraints are present, this symmetry implies the conservation of linear 
momentum. In the presence of nonholonomic constraints, a modified conservation law 
holds: it is easy to see that the following vector field annihilates F along C: 

ox oy 09 

(Any scalar multiple of the above vector field is also allowed.) This generalized vector 
field corresponds with the section £ = (—Ry',Rx', 1) of tt\ Q g F . As £ is vertical, the 
nonholonomic momentum lemma [4~9l can be applied. 

The right-hand side of (H]) is 

j 2 <t>*^(Lr)) = [ - Rpy'x + Rpx'y - KRx"'y" + KRy"'x"]n. 
The nonholonomic momentum map J ', on the other hand, is given by 

Jf h ' = - [p{Rx'y - Ry'x) + a9]ds - [KR(y'x"' - x'y'") + 0]dt, 
and the nonholonomic momentum equation hence becomes 

Ry'(px + Kx"") - Rx'{py + Ky"") = a9 - (36". (21) 



This conservation law can also be derived from the nonholonomic field equations (|17|) 
by subtracting the second equation multiplied by x' from the first equation multiplied 
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by y' , and using the third equation to eliminate the Lagrange multipliers A and [i. 
Unfortunately, the knowledge of this nonholonomic conservation law does not help us 
in solving the field equations (in contrast with the situation for the vertical rolling 
disc; see [4]). 

6. Conclusions 

In this paper, we presented a geometric framework for nonholonomic field theories 
with symmetries. We showed that the momentum map associated with a group 
action satisfies a certain momentum equation, which we proved in a number of cases. 
On the one hand, there exists a momentum equation for nonholonomic symmetries 
which act nontrivially on the base space (examples being energy conservation for the 
nonholonomic rod), while on the other hand a similar result exists for generalized 
symmetries associated to a vertical group action. The Bcnenti system is an example 
of the latter. 

It is likely that the results in this paper can be generalized still further. In particular, 
there seems to be no reason why there shouldn't be a momentum equation for 
generalized symmetries associated to a non- vertical group action, up to some technical 
restrictions (for instance, that the generalized vector field should be related to a regular 
vector field on the base space) . This generalization would encompass both momentum 
equations derived in this paper. Moreover, it would be interesting to study further 
examples of such a momentum equation. 
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Appendix: noncovariant nonholonomic constraints 

A special class of nonholonomic field theories consists of those where the base space 
X is R x M, where the first factor represents time, and such that ir is trivial. In other 
words, there exists a canonical distinction between time and space. Accordingly, one 
can show that in this case, the jet bundle J 1 ^ is isomorphic to the product bundle 
K x [J 1 (M, S) x 5 TS] , thus providing a canonical distinction between derivatives of 
the fields with respect to time and space. 

This class of field theories was discussed in detail in [23] and includes among others also 
theories of nonrelativistic elasticity (see also [17]). The example studied in section I5T21 
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is also among these field theories. 

As shown in [23], the reaction forces for this kind of field theory are not the ones 
obtained from the Chctaev prescription in remark T3.6l The problem is that this form 
of the principle is fully covariant, in the sense that no distinction is made between 
spatial derivatives and derivatives with respect to time. However, in some cases 
(the nonholonomic Cosserat rod being one of them) the time derivatives do play 
a distinguished role, and therefore a different, "noncovariant" Chetacv principle is 
needed. 

The main objective of this appendix is to propose such a principle, using the geometry 
of J 1 ^, and especially the isomorphism with M. x [J 1 (M, S) x sTS] . Using this modified 
principle, we derive the correct form of the reaction forces for the nonholonomic 
Cosserat rod, which were obtained in [23] using a different approach. 



1.1. A new vertical endomorphism 

Recall the coordinate expression |T]) of the vertical endomorphism S v on J lr K. This 
tensor field was constructed by Saunders [20] using a map assigning to each one-form 
uj on X the vector-valued one-form on J 1 ^ given in coordinates by (see [20, p. 156]) 

= u fi (dy a -y^x")®-^, where u^c^daA (22) 

Roughly speaking, the vertical endomorphism then arises, once a volume form on 
X is chosen, by putting 

S„ = Sd^A^d"^), (23) 

where the wedge operator 'A' is defined as follows: if $ is a vector- valued fc-form 
on J 1 ^, and a is a regular (i.e. R- valued) Worm, then $Aq is the vector- valued 
(k + Z)-fbrm given by ($Aa, f3) = ($, (3) A a for all (3 e f2 1 (J 1 7r). 

It is obvious that (|23|) is fully covariant, in the sense that no distinction is made 
between the variables on the base space. However, in elastodynamics, this is not 
always desirable, as we have a distinguished direction of time. Therefore, we propose 
the following "non-covariant" vertical endomorphism: 

Definition 1.1. The non-covariant vertical endomorphism is the vector-valued (n+1)- 
form S n . c . defined as S n . c . := Sdt ^{^i'Hm), where Sdt is the vector-valued one-form 
associated to dt as in i22\) . 

Here t\m denotes the volume form on M defined by restriction to M of d n Xo 

Note that dt is a well defined one-form on X = M x M; therefore, S n . c . is an intrinsic 
object. In coordinates, S n . c , is given by 

Sn.c. = (dy a - i/«dx") A d"a; ® A 

Remark 1.2. In [20], it is shown that the action of the vertical endomorphism S v 
is related to the fact that tti^q : J 1 ir — > Y is an affine bundle; i.e. there exists an 
affine action of ir*T*X ® Vir on J l ir. The noncovariant vertical endomorphism can 
be understood in a similar vein, by restricting this affine action to ir*T*M. <g> Vtt. o 
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1.2. The bundle of constraint forces 

Let i : C c — > J X 7r be a constraint manifold. In section [31 we required that reaction 
forces be n-horizontal and 1-contact. This leads to local expressions of the form ([4J. 
For noncovariant constraints, we attribute a special status to the time coordinate, and 
therefore we require also that the following holds: 

i v i w $ = 

for all tangent vectors v,w on J 1 ^ such that T(pr 1 o tti)(v) = T(pr x o ni)(w) = 0, 
where pr x : R x M — > R is the projection onto the first factor (this condition expresses 
that v and w do not contain a component proportional to Jj). In coordinates, this 
implies that <£> has the following form: 

<S> = A a (dy a -y«dx»)Ad n x , (24) 

where the A a are local functions on C. Compare this with (j4|) and note that only the 
"time" component is left [i.e. the component proportional to 6 a A (J^J r/)). 

Using the noncovariant vertical endomorphism of section fl. 11 one can construct from 
C a natural candidate for the bundle F. This is a field-theoretic generalization of the 
well-known Chetaev principle from mechanics. 

Recall that C is assumed to be given by the vanishing of k functions ip a and define the 
associated bundle of reaction forces as the bundle F locally spanned by the following 
(n + l)-forms: <!>" := S* c (dip a ), or in coordinates: 

$« = |g(dy a -^d^)Ad"zo. (25) 

The (n + l)-form is therefore of the form outlined in ([23]). with A" = f^. 
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